Multilevel Variable-Block Schur-Complement-Based Preconditioning for the Implicit Solution of the Reynolds- Averaged Navier-Stokes Equations Using Unstructured Grids

Implicit methods based on the Newton’s rootfinding algorithm are receiving an increasing attention for the solution of complex Computational Fluid Dynamics (CFD) applications due to their potential to converge in a very small number of iterations. This approach requires fast convergence acceleration techniques in order to compete with other conventional solvers, such as those based on artificial dissipation or upwind schemes, in terms of CPU time. In this chapter, we describe a multilevel variable-block Schur-complement-based preconditioning for the implicit solution of the Reynolds-averaged Navier-Stokes equations using unstructured grids on distributed-memory parallel computers. The proposed solver detects automatically exact or approximate dense structures in the linear system arising from the discretization, and exploits this information to enhance the robustness and improve the scalability of the block factorization. A complete study of the numerical and parallel performance of the solver is presented for the analysis of turbulent Navier-Stokes equations on a suite of three-dimensional test cases

preliminary design of a radial turbine for methane expander rocket engine

Abstract The present paper summarizes recent research efforts carried out at the Italian Aerospace Research Centre, CIRA, aimed at using radial turbines in modern rocket-propulsion systems. Over the last few years, CIRA has been involved in the HYPROB program, funded by the Italian Ministry of Research (MIUR). HYPROB aims at developing competences to consolidate the national background on rocket-engine systems for future space applications. Since nowadays liquid methane represents an innovative fuel in aerospace propulsion, one of HYPROB's objectives is the development of a simulation tool for the preliminary design of the radial turbines used to drive the turbo-pumps in expander-cycle rocket-engines operating with methane

Generation of Surface Maps of Erosion Resistance for Wind Turbine Blades under Rain Flows

Rain erosion on wind turbine blades raises considerable interest in wind energy industry and research, and the definition of accurate erosion prediction systems can facilitate a rapid development of solutions for blade protection. We propose here the application of a novel methodology able to integrate a multibody aeroelastic simulation of the whole wind turbine, based on engineering models, with high-fidelity simulations of aerodynamics and particle transport and with semi-empirical models for the prediction of the damage incubation time. This methodology is applied to generate a parametric map of the blade regions potentially affected by erosion in terms of the fatigue life of the coating surface. This map can represent an important reference for the evaluation of the sustainability of maintenance, control and mitigation interventions

Preliminary Design Method of a Turbopump Feed System for Liquid Rocket Engine Expander Cycle

Abstract The present research effort deals with simplified theoretical models for the preliminary design and performances assessment of centrifugal pumps for liquid rocket propulsion. These models have been developed within the Concurrent Design Facility, under development at the Italian Aerospace Research Centre (CIRA), in the framework of the HYPROB program. In particular, this work is aimed at developing a theoretical model, via the implementation of a MatLab code, capable to predict the geometry and performance of centrifugal turbopumps, thus providing useful indications for the preliminary design of the turbopump feed system

The role of mesh generation, adaptation, and refinement on the computation of flows featuring strong shocks

Within a continuum framework, flows featuring shock waves can be modelled by means of either shock capturing or shock fitting. Shock-capturing codes are algorithmically simple, but are plagued by a number of numerical troubles, particularly evident when shocks are strong and the grids unstructured. On the other hand, shock-fitting algorithms on structured grids allow to accurately compute solutions on coarse meshes, but tend to be algorithmically complex. We show how recent advances in computational mesh generation allow to relieve some of the difficulties encountered by shock capturing and contribute towards making shock fitting on unstructured meshes a versatile technique

A new computational technique for re-entry flow calculations based upon a shock-fitting technique for unstructured grids

An in-house developed, 2D/3D unstructured CFD solver has been extended to deal with a mixture of thermally perfect gases in chemical non-equilibrium. The Euler equations have been coupled with a state-to-state kinetic model for argon plasma. The spatial discretization uses compact stencil Residual Distribution Schemes and shock waves can be modelled using either shock-capturing or shock-fitting. Promising results have been obtained using the shock-fitting approach for a 2D hypersonic flow past the fore-body of a circular cylinder

Assessing wind turbine energy losses due to blade leading edge erosion cavities with parametric CAD and 3D CFD

Wind turbine leading edge erosion is a complex installation site-dependent process that spoils the aerodynamic performance of wind turbine rotors. This gradual damage process often starts with the formation of pits and gouges leading ultimately to skin delamination. This study demonstrates the application of open source parametric CAD functionalities for the generation of blade geometries with leading edge erosion damage consisting of pits and gouges. This capability is key to the development of high-fidelity computational aerodynamics frameworks for both advancing knowledge on eroded blade aerodynamics, and quantifying energy losses due to erosion. The considered test case is an offshore 5 MW turbine featuring leading edge pit and gouge damage in the outboard part of its blades. The power and loads of the nominal and damaged turbines are determined by means of a blade element momentum theory code using airfoil force data obtained with 3D Navier-Stokes computational fluid dynamics. An annual energy loss between about 1 and 2.5 percent of the nominal annual energy yield is encountered for the considered leading edge damages. The benefits of adaptive power control strategies for mitigating erosion-induced energy losses are also highlighted

Shock-fitting and predictor-corrector explicit ALE Residual Distribution

International audienc

Using the VBARMS method in parallel computing

The paper describes an improved parallel MPI-based implementation of VBARMS, a variable block variant of the pARMS preconditioner proposed by Li, Saad and Sosonkina [NLAA, 2003] for solving general nonsymmetric linear systems. The parallel VBARMS solver can detect automatically exact or approximate dense structures in the linear system, and exploits this information to achieve improved reliability and increased throughput during the factorization. A novel graph compression algorithm is discussed that finds these approximate dense blocks structures and requires only one simple to use parameter. A complete study of the numerical and parallel performance of parallel VBARMS is presented for the analysis of large turbulent Navier-Stokes equations on a suite of three- dimensional test cases